AD-A079  962 
UNCLASSIFIED 


KENTUCKY  UNIV  LEXINGTON  OE*»T  OF  ELECTRICAL  ENGINEERING  F/G  20/14 
ON  THE  PHYSICAL  REALIZA3ILITY  OF  BROADBAND  LUMPED-PARAMETER  EOU— ETC(U) 
JUL  79  L  *  PEARSON*  D  R  HILTON  AF0SR-77-3452 

RR-79-2  AF0SR-TR-B0-0004  *  NL 


zammszmxm 


PORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


l.  GOVT  ACCESSION  NO.I  1.  RECIPIENT’S  CATALOG  NUMBER 


THE^HYSICAL  REALIZABILITY  OF  BROADBAND 
^  .>JDMPEI>^ARAMETER  EQUIVALENT  jgLRCUITS  FOR  / 


S.  TYPE  OP  REPORT  a  PERIOD  COVERED 


ENERGY-' 


COLLECTING  STRUCTURES, 


7.  AUTHOR^ 


O^O.  REP^T  NUMBER 


9.  CONTRACT  OR  GRANT  HUMBERT*; 


iLSsrfe,7  (Sl  ^I'E  Q?}-  ^AMki77-3452 


msjMm i 


*.  perporming  organization  name  and  address 

Department  of  Electrical  Engineering  ' 
University  of  Kentucky 
Lexington.  KY  40506 


<1.  CONTROLLING  OPPICE  NAME  ANO  ADDRESS 

Physics  Directorate 

Air  Force  Office  of  Scientific  Research 
Bolling  AFB,  Washington,  DC  20332 


«.  MONITORING  AGENCY  NAME  »  AOORESSflf  dillmront  Inm  Conlnllii-d  Offfo)  IS.  SECURITY  CLASS,  (of  Ml.  report) 

qjrQdJ  unclassified 

IS*.  OECLASSIPICATION/ DOWNGRADING 
SCHEDULE 


tc.  DISTRIBUTION  STATEMENT  (of  Ml.  Ropo.l) 


MBER  OF  PAGES 

41 


Approved  for  public  release j 
distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (of  Mo  aAolroct  antorod  in  Block  20,  1 1  dIMoront  from  Report) 


It.  KEY  WORDS  (Conttnum  on  revtnt  II  noooooorr  md  Identify  by  Woe*  numb*) 

Transient,  equivalent  circuits,  antenna,  scattering,  coupling, 
singularity  expansion 


20.  ABSTRACT  'Continue  on  rtvorao  eldo  It  noefieorr  and  Identity  by  block  number) 

*  The  practical  implementation  of  equivalent  circuits  which,  represent  the  energy 
-  delivered  to  an  arbitrary  load  Impedance  by  an  antenna  or  a  port  on  a  scatterer 
is  discussed.  Particular  attention  is  paid  to  the  positive-real  function 
issues  assdeiated  with  the  equivalent  admittance  of  the  structure  when  evolved 
from  the  point  of  view  of  the  Singularity  Expansion  Method.  A  proof  is  given 
which  establishes  that  the  admittances  associated  with  current  eigenmode  con¬ 
tributions  to  the  total  current  are  individually  positive-real  functions.  An  - 


I  WLW.IeMA 


ft 


20.  Abstract  continued 

ad  hoc  approach  to  approximating  these  admittances  with  simpler 
pole-pair  positive-real  admittances  is  given  for  high-Q  structures. 
A  brief  survey  of  the  applicable  circuit  synthesis  algorithms  and 
their  features  is  provided.  Some  fundamental  considerations  in 
the  implementation  of  the  sources  representing  the  energy  coupled 
from  the  incident  wave  are  presented.  ^ 


1 


1 


ON  THE  PHYSICAL  REALIZABILITY 
OF  BROADBAND  LUMPED-PARAMETER 
EQUIVALENT  CIRCUITS  FOR 
ENERGY-COLLECTING  STRUCTURES 


L.  Wilson  Pearson 

Department  of  Electrical  Engineering 
University  of  Kentucky 
Lexington,  KY  40506 


Donald  R.  Wilton 

Department  of  Electrical  Engineering 
University  of  Mississippi 
University,  MS  38677 


University  of  Kentucky 
Electromagnetics  Research  Report 
79-2  / 

July,  1979 


D 


Fiscal  Report  AFOSR  Grant  No. 


AFOSR  77-3452. 

‘-S.'vARCH  (AFSC) 


appro  .•  a 

Diilri  i0  a-i 

A.  D • 


Technical  n,.,-  .--.tlon  OiYicur 


-  •‘'•■d  and  ig 

(7b). 


1 


1  f? 


ABSTRACT 


ON  THE  PHYSICAL  REALIZABILITY  OF  BROADBAND  LUMPED-PARAMETER 
EQUIVALENT  CIRCUITS  FOR  ENERGY-COLLECTING  STRUCTURES 

The  practical  Implementation  of  equivalent  circuits  which  represent  the 
energy  delivered  to  an  arbitrary  load  impedance  by  an  antenna  or  a  port  on  a 
scatterer  is  discussed.  Particular  attention  is  paid  to  the  positive-real 
function  issues  associated  with  the  equivalent  admittance  of  the  structure 
when  evolved  from  the  point  of  view  of  the  Singularity  Expansion  Method.  A 
proof  is  given  which  establishes  that  the  admittances  associated  with  cur¬ 
rent  eigenmode  contributions  to  the  total  current  are  individually  positive- 
real  functions.  An  ad  hoc  approach  to  approximating  these  admittances  with 
simpler  pole-pair  positive-real  admittances  is  given  for  hlgh-Q  structures. 

A  brief  survey  of  the  applicable  circuit  synthesis  algorithms  and  their 
features  is  provided.  Some  fundamental  considerations  in  the  implementation 
of  the  sources  representing  the  energy  coupled  from  the  incident  wave  are 
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ON  THE  PHYSICAL  REALIZABILITY  OF  BROADBAND  LUMPED-PARAMETER 
EQUIVALENT  CIRCUITS  FOR  ENERGY-COLLECTING  STRUCTURES 

I.  INTRODUCTION 

Schelkunoff  presented,  to  our  knowledge,  the  first  development  relative 
to  physically  realizable  equivalent  circuits  for  radiating  structures  [1]. 
While  his  theory  offers  a  great  deal  of  Insight  into  the  resonance  and  anti- 
resonance  phenomena  which  antennas  manifest,  the  theory  provides  no  means  for 
the  actual  physical  determination  of  equivalent  circuits  for  radiating  or 
receiving  structures.  More  recently,  Weber  and  Toulios  provided  a  means  for 
some  equivalent  circuit  modeling  of  antennas  on  a  broadband  transient  basis 
by  using  a  Laplace  domain  curve-fitting  technique  [2]. 

In  1971,  Baum  introduced  the  Singularity  Expansion  Method  (SEM)  [3]. 

SEM  provides  a  natural  resonance  representation  for  the  current  on  an  antenna 
or  scattering  structure.  Such  a  complex  resonance  representation  admits  to 
the  possibility  of  applying  general  formal  circuit  synthesis  procedures  which 
can  be  implemented  in  specific  cases,  provided  the  SEM  description  for  the 
case  in  question  is  available.  However,  a  number  of  questions  have  arisen 
relative  to  the  completeness  of  the  representation.  A  great  deal  of  work  in 
the  SEM  area  has  proceeded  since  Baum's  introduction  of  it.  The  three  sum¬ 
mary  works  by  Baum  [4,5,6]  provide  excellent  bibliographies  and  give  a  nearly 
up-to-date  summary  of  what  has  proceeded.  In  a  later  report  [7],  Baum 
relates  the  singularity  expansion  to  eigenfunction  analysis  of  integral 
operators.  This  relationship  bears  heavily  on  the  work  reported  here.  Also 
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in  a  later  note,  he  develops  a  formal  theory  for  equivalent  circuit  synthesis 
for  radiating  and  energy  collecting  structures  [8]. 

Hess  has  used  SEM  analysis  to  develop  a  transfer  function  model  for 
determining  EMP  coupling  by  means  of  equivalent  transfer  function  circuits 
[9].  The  present  method  goes  beyond  the  work  of  Hess  in  that  it  accounts 
for  the  possibility  of  loading  at  the  terminals  of  the  receiving  structure. 
More  recently,  Schaubert  has  used  measured  transient  response  data  which  was 
subsequently  analyzed,  using  a  Prony-type  algorithm,  for  the  determination 
of  the  complex  natural  resonances  in  a  network  synthesis  [10].  This  latter 
work,  though  ad  hoc  in  character,  is  likely  to  be  quite  important  in  elabora¬ 
tions  on  the  synthesis  methods  presented  here.  The  use  of  measured  data 
allows  the  possibility  of  structures  which  are  too  complex  to  be  modeled  on 
a  theoretical  or  numerical  basis.  Schaubert,  however,  limits  his  concern  to 
the  determination  of  the  admittance  or  impedance  part  of  the  equivalent  cir¬ 
cuit;  he  does  not  consider  the  source  problem. 

In  this  work,  we  have  explored  some  of  the  ramifications  of  Baum's 
development  in  reference  [8].  In  particular,  we  have  emphasized  the  physi¬ 
cal  realizability  of  the  circuits  developed  with  passive  elements.  As  a 
consequence,  the  dominant  issue  which  arises  is  that  of  the  positive-realness 
(PR-ness)  of  the  impedance  or  admittance  representations  which  result.  In  a 
companion  work  [11],  the  results  of  an  extensive  numerical  study  emphasizing 
PR  properties  of  admittances  for  straight  wires  and  wire  loop  structures  are 
presented.  In  subsequent  paragraphs  here,  we  shall  draw  on  the  results  of 
reference  [11]  from  time  to  time  as  evidence  supporting  various  assertions. 

In  the  following,  we  develop  the  formal  theory  of  equivalent  network 
synthesis  based  on  the  singularity  expansion.  This  development  essentially 
follows  Baum  [8].  It  is  reiterated  here  for  the  sake  of  completeness  and 
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for  the  sake  of  introducing  what  we  feel  to  be  a  somewhat  simplified  nota¬ 
tion.  We  have  chosen  to  limit  our  discussion  to  that  of  synthesis  of  the 
admittance  of  the  antenna  or  energy  collector,  and  to  electric  field  inte¬ 
gral  equation  formulations.  These  results  may  be  generalized  to  embrace 
impedance-type  synthesis  and  also  alternative  integral  equations  for  the 
structures.  Indeed,  Baum's  formal  development  in  reference  [8]  Includes 
this  level  of  generality.  We  have  chosen,  for  the  sake  of  clarity,  t:o 
restrict  the  scope  of  the  present  development. 


II.  FORMAL  DEVELOPMENT  OF  THE  EQUIVALENT  CIRCUIT:  TERMINAL  EIGENELEMENTS 

In  this  section  we  develop  a  formal  equivalent  circuit  for  an  energy¬ 
collecting  structure.  By  an  energy-collecting  structure,  we  mean  either  a 
single-port  antenna  or  a  scatterer  with  an  inadvertent  penetration  which  can 
be  identified  in  terms  of  a  port. 

The  objective  of  the  development  is  indicated  in  Figure  la.  There,  a 
generic  electromagnetic  energy  collecting  structure,  loaded  by  a  frequency 
dependent  admittance  Y^(s),  is  indicated.  (We  shall  use  the  complex 
frequency  variable  s— the  Laplace  transform  variable— to  indicate  frequency 
thoughout  this  development.  A  tilde  ~  over  a  function  denotes  a  Laplace 
transformed  quantity.)  We  wish  to  replace  the  influence  of  the  driven 
energy  collecting  structure  on  the  admittance  Y^(s)  with  a  Norton  equivalent 
circuit,  as  shown  in  the  right  portion  of  Figure  la.  The  elements  in  the 
equivalent  circuit  are  derivable  by  the  steps  indicated  in  Figures  lb  and 
lc.  First,  the  equivalent  admittance  of  the  structure  Y  (s)  is  obtainable 

1j 

by  treating  the  structure  as  a  transmitting  antenna.  To  do  this,  we  intro¬ 
duce  a  frequency  dependent  voltage  V(s)  across  the  port  of  the  structure, 
determine  the  port  current  I(s),  and  find  the  ratio  I(s)/V(s).  This  ratio 
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Figure  1.  Development  of  the  Norton  equivalent  circuit  for  an  energy 
collecting  structure:  (a)  the  equivalence  sought;  (b)  the  admittance 
problem;  and  (c)  the  current  source  problem. 
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is  the  admittance  in  question.  The  source  term  in  the  equivalent  circuit 

I 

is  determined  by  shorting  the  port  and  treating  the  resulting  scattering 

problem.  From  it,  we  obtain  the  net  current  at  the  port,  I  (s) ,  which 

sc 

becomes  the  generator  in  the.  equivalent  circuit. 

The  formal  SEM  equivalent  circuit  results  from,  representing  the  solu¬ 
tions  to  the  two  boundary  value  problems  indicated  in  Figures  lb  and  lc  in 
terms  of  the  singularity  expansion.  In  the  present  development,  we  use  the 
eigenfunction  expansion  for  the  solutions  to  the  respective  problems  as  an 
intermediate  step  coward  the  development  of  the  singularity  expansion.  The 
intermediate  forms  of  the  eigenfunction  expansions  and  properties  associated 
therewith  subsequently  prove  useful  in  the  development  of  practical  equiva¬ 
lent  circuits. 

Either  of  the  boundary  value  problems  in  question  may  be  formally  cast 
as  an  integral  equation  for  the  surface  current  on  the  surface  of  the  energy 
collecting  structure.  The  electric  field  integral  equation  for  the  struc¬ 
ture  is  representable  as 

<f(r,r',s)  ;  J(r',s)>  -  Einc(r,s)  ,  r  e  S  ,  (1) 

where  the  dyadic  kernel  is  given  by 

r(r,r'  ,s)  -  [spl  -  ~  grad  div]  e”«  I  r~r'  I /c/  [Arrl  r-r*  |  /c]  ,  (2) 

SC 

with 

u  ■  permeability  of  the  medium, 
e  ■  permittivity  of  the  medium,  and 

_  i. 

c  ■  (pc)  2  -  velocity  of  light  in  the  medium. 
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-inc  - 

The  tern  E  (r,s)  represents  the  tangential  component  of  the  incident  field 
on  s  exciting  the  object.  The  bracket  notation  is  used  to  indicate  a  sym¬ 
metric  product— a  surface  integral  over  the  extent  of  the  scatterer  surface, 
viz. , 

<r(r,r' ,s)  ;  J(r’ ,s)>  «  //  r(r,r',s)  *  J(r*,s)  dS'  .  (3) 

S 

The  term  E  (r,s)  is  replaced  with  a  field  due  to  the  impressed  voltage 
V(s)  in  the  case  of  the  transmitting  problem  of  Figure  lb. 

A  formal  solution  to  the  integral  equation  (1)  >>ove  is  given  in  terms 
of  the  eigenfunction  expansion  for  the  surface  current  J(r,s): 


J(r’3)  “  l  Tm 


<J  (r,s)  ;  E  c(r,s)> 

— r — - — r — -  J  (r,s)  . 

<Jn(r,s)  ;  Jn(r,s)> 


(4) 


The  Jn(r,s)  and  the  XQ(s)  are  the  eigenfunction/eigenvalue  pairs  for  the 
integral  operator.  Namely, 


<r(r,r',s)  ;  Jn(r’,s)>  -  Xn(s)  Jn(r,s)  . 


(5) 


The  admittance  of  the  structure  is  determined  by  solving  the  antenna 
problem  indicated  in  Figure  lb.  The  incident  field  is  taken  to  be  the 
quasi-static  EMF  across  the  gap  width  A  due  to  the  impressed  voltage  V(s), 

ginC(  a)  s  v(s) / A  a  ,  r  e  S  , 

8  g 

with  a  a  unit  vector  oriented  from  one  gap  surface  to  the  other  on  S  , 
g  6  g* 

defining  the  reference  direction  for  the  EMF,  as  depicted  in  Figure  2.  In 
terms  of  the  eigenfunction  expansion,  this  solution  may  be  written  as 
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Jant(r,s) 


lisl  z  _JL_ 

A  n  X  (a) 


<\G'.s)  ;  ag>g 


n 


<Jn(r,s)  ;  Jn(r,a)> 


-  .  *■  J  <r,s)  • 


(6) 


The  S  on  Che  bracket  of  Che  coupling  integration  in  the  numerator  indicates 
8 

that  the  extent  of  the  integration  is  limited  to  the  surface  area  comprising 
the  gap  region  of  the  structure.  The  net  current  flowing  through  the  gap 
region  is  representable  as  an  azimuthal  integration  around  the  gap  of  the 
expression  (6).  This  integration  may  be  stated  formally  as  1/A  times  the 
integration  over  Sg.  Thus, 

I(s)  =  j  <ag  ;  3ant(r,s)>s 

s 


r 

A2  n 


,  4n(^s)  >  Vsg 

*n^  <Jn(r,s)  ; 


The  eigenfunction  expansion  of  the  admittance  Y(s)  follows  directly: 


Y(s)  -  I(s)/V(s) 


■  l 
n 


<Jn(r,s)  ;  ag>s 


£_ 


Xn^  A2  <Jn(r,s)  ;  JQ(r,s)> 


(7) 


£  Yn<S) 

n  n 


We  introduce  the  notation  Y^Cs)  to  represent  the  elements  of  the  summation 
in  equation  (7).  We  term  these  Yn  "terminal  eigenadmittances".  We  use  this 
terminology  to  explicitly  distinguish  these  terms  from  the  eigenadmittance 
quantities  which  Baum  defines  in  reference  [7]. 


The  short  circuit  current  for  the  boundary  value  problem  of  Figure  1c 
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The  expressions  above  lead  to  the  equivalent  circuit  representation 

shown  in  Figure  3.  This  circuit  is  a  recasting  of  the  equivalent  circuit 

sought  in  Figure  la.  The  recasting  results  from  pairing  the  current  sources 

and  the  admittance  contributions  due  to  individual  elements  in  the  eigenset. 

The  form  of  the  circuit  suggests  transforming  from  a  Norton  to  a  Thevenin 

form  for  each  terminal  eigensource/eigenadmittance  pair,  as  shown  in  Figure 

3b.  The  elements  of  this  circuit  are  given,  respectively,  by  equations  (7) 

and  (11) .  A  caution  must  be  given  in  regard  to  applying  these  forms.  Note 

that  if  the  symmetric  product  <Jn  ;  a^>  *  0*  Y Q  vanishes  quadratically  in 

**  8 

this  product  while  has  a  reciprocal  singularity.  Thus,  for  such  a  ter¬ 
minal  location,  although  Vn  -*■  «,  the  entire  leg  of  the  circuit  vanishes. 


III.  FORMAL  DEVELOPMENT  OF  THE  EQUIVALENT  CIRCUIT:  SINGULARITY  EXPANSION 
OF  ELEMENTS 

The  singularity  expansion  for  surface  current  quantities  on  a  scatter¬ 
ing  object  can  be  obtained  from  the  eigenfunction  expansion  by  performing 
the  residue  series  expansion  of  expressions  such  as  equation  (4)  in  terms 
of  the  poles  or  complex  natural  resonances  of  the  object.  These  poles  are 
the  complex  frequencies  at  which  the  eigenvalues  of  the  integral  operator 
vanish,  viz.. 


(3ni  !  Vsni> 


0}  . 


Recently,  we  have  uncovered  evidence  that  branch  points  may  be  present  in  the 
complex  s-plane,  and  that  at  least  two  eigenvalues  are  degenerate  at  these 

frequencies.  Marin  and  Latham  [12)  argue  the  absence  of  branch  points  in  the 

« 

complex  plane  representation  of  the  eigenvalues.  However,  their  argument  is 
based  on  an  initial  assumption  of  distinct  eigenvalues,  whereas,  in  fact,  the 


equivalence  circuit  organized  into  eigenelement  pairs 
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branch  points  actually  originate  from  the  eigenvalue  degeneracies  [13]. 
However,  they  correctly  argue  that  for  an  object  Inmersed  In  a  lossless 
medium  the  overall  representation  of  the  current  Is  free  of  branch  Integral 
contributions.  There  appears  to  be  no  Inconsistency  here,  for  in  the  sum¬ 
mation  over  all  the  eigenvalues  as  Indicated  In  equation  (4) ,  the  branch 
integral  contributions  contributed  by  all  of  the  eigenvalues  associated  with 
a  given  branch  sum  to  zero.  The  evidence  of  such  a  phenomenon  occurring  In 
singularity  expansion  representations  for  scattering  objects  is,  at  this 
writing,  sketchy,  although  we  plan  to  elaborate  on  this  subject  in  the 
future.  Because  the  presence  of  branch  integrals  in  individual  eigenvalue 
representations  is  potentially  highly  significant  in  terms  of  the  physical 
realizability  of  the  admittances  shown  in  Figure  3,  we  shall  explicitly 
include  branch  integral  contributions  in  our  admittance  representations  for 
the  moment. 

The  reciprocal  eigenvalue  factor  which  appears  in  equation  (4)  may  be 
expanded  in  a  residue  series  expansion  [6]  to  yield 


An(s) 


+  b  (s)  +  e  (s)  . 
n  n 


(12) 


The  residues  in  this  expansion  are  given  by 
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This  form  results  from  applying  residue  calculus  to  the  Rayleigh  quotient 
representation  of  the  eigenvalues.  The  bn(s)  are  contributions  to  the 

me 

eigenvalue  expansion  arising  from  branch  integrals  and  the  en(s)  are  entire 
functions  which,  according  to  the  Mittag-Leffler  theorem,  may  also  be 
required.  Neither  of  these  types  of  functions  are  represented  explicitly 
here,  but  are  included  in  the  expansion  above  for  the  sake  of  completeness. 
It  appears,  to  date,  that  for  practical  circuit  synthesis  it  is  necessary 
that  both  of  these  terms  be  zero.  No  completely  conclusive  arguments  have 
been  put  forth  that  they  are  zero  in  the  general  situation. 

The  use  of  the  expansion  (12)  and  the  residue  definition  (13)  in  the 
terminal  eigenadmittances  of  (7)  yields  the  following: 


Y 

n 


(s)  =  Z  B 


i 


ni 


<J  (r,s  .)  ;  a 
“  ni'  g  sr 


1 


s  -  s 


ni 


+ 


+  ?i(s)  +  Y*nt(s) 
n  n 


(14) 


where 


6ni  “  <3n(9*sni)  ;  ’8ni)/ds  *  3n(?'»,ni)>‘1 


The  Yq  and  the  Y^  are,  respectively,  admittance  contributions  from  the 
branch  integral  function  bQ  in  (12)  and  from  the  entire  function  eQ.  We 
emphasize  that  for  numerically  or  experimentally  derived  SEM  descriptions 
neither  of  these  terms  are  explicitly  Identifiable.  By  defining 


Ani  “  6ni 


<J  (r,s  .)  ;  a 
n  »i  8  Sr 


we  may  compactly  write  the  expansion  (14)  as 


V8> 


l  Yni<8> 


(15) 


The  branch  integral  and  entire  function  terms  are  dropped  here  and  henceforth 
in  expressions  for  Yq.  The  potential  consequences  of  the  presence  of  a  non¬ 
zero  Yt)i  are  discussed  in  Section  IV.  The  Y  .  (s)  are  termed  modified  pole 
n  ni 

admittances*.  We  note  that  J  , (r.s  ,)  solves 

ni  *  ni 

<*(r»r’ ,snl)  ;  3n(r',snl)>  -  0 

and  hence  is  the  natural  mode  of  the  object  associated  with  the  pole  s  . . 

m 

It  is  evident  from  (10)  and  (11)  that  the  terminal  eigensource  currents 
may  be  expanded  in  the  same  fashion  to  yield 

I8c(s)  »EV(s.)A,  T - - -  +  —1 

n  i  n  ni  ni  [s  -  sq1  «nlJ 

+  1^(3)  +  I®nt(s) 
n  n 

-  I  \i  *  lnl<»  +  C'O)  •  «« 

We  term  the  Vni  pole  source  voltages.  Through  precisely  the  same  reasoning 
as  that  used  on  the  eigensource  current/eigenadmittance  pairs,  we  interpret 
(15)  and  (16)  as  representing  parallel  combinations  of  legs  comprising 

«•  «M 

series  connections  of  Y  , (s)  and  V  J(s). 

ni  ni 


*We  use  the  term  "modified"  to  be  consistent  with  Baum  in  [8).  There  he 
defines  an  "unmodified"  pole  admittance  A./(s-s.),as  well. 

I1 1  * 
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IV.  POSITIVE-REAL  CONSIDERATIONS 

The  fundamental  condition  which  dictates  the  physical  realizability  of 
any  impedance  or  admittance  is  that  it  be  a  positive-real  (PR)  function  of 
the  complex  frequency.  This  condition  results  from  the  fact  that  a  passive 
structure  (or  network)  can  never  dissipate  a  negative  amount  of  energy. 

Since  the  admittance  function  (7)  represents  that  manifested  by  a  passive 
energy-collecting  object,  it  must  be  a  PR  function.  If  this  admittance  is 
to  be  approximately  realized  by  a  lumped-parameter  circuit,  it  is  necessary 
that  the  summation  of  (7)  be  truncated  to  some  finite  n.  Further,  each 
Yq(s)  embraces  a  summation  over  a  potentially  infinite  set  of  poles,  as 
expressed  in  (14).  This  summation,  too,  must  be  truncated  for  circuit 
realization.  The  circuit  realization  will  encompass  a  combination  of  cir¬ 
cuits  representing  terms  in  either  the  summation  of  (14)  or  of  (7).  There¬ 
fore,  two  issues  emerge  relative  to  the  way  in  which  this  approximate 
realization  is  approached:  1)  Is  it  possible  to  truncate  the  summation (s) 

of  (14)  and  (7)  in  such  a  way  as  to  maintain  PR  properties,  and  2)  are  the 

«* 

elemental  terms  in  the  suranand  of  (14)  or  the  Y  (s)  of  (7)  individually  PR? 

n 

This  section  deals  with  these  issues,  as  we  presently  understand  them. 

The  terminal  elgenadmittances  Y^  (s)  are  positive-real  functions  of 
frequency.  This  result  is  developed  in  detail  in  the  Appendix.  The  develop¬ 
ment  given  there  follows,  in  essence,  that  which  one  uses  in  developing  PR 
properties  for  circuits.  The  additional  feature  which  must  be  brought  to 
bear  in  the  distributed  problem  of  an  energy  collector  is  the  spatially- 
selective  feature  of  the  eigenmode  currents. 

That  the  terminal  elgenadmittances  are  positive-real  functions  of  fre¬ 
quency  is  a  particularly  robust  result.  First  of  all,  it  provides  a  means 
for  truncating  the  summation  given  in  (7).  Since  any  terminal  eigenadmittance 
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term  Is  by  Itself  PR,  then  any  finite  sum  of  terminal  eigenadmittances  Is 
also  PR.  Thus,  the  summation  in  (7)  may  be  truncated  to  produce  a  PR 
approximation  to  Y(s).  The  truncation  may  be  based  either  on  temporal  or 
on  spatial  frequency  of  the  excitation.  The  lumped-parameter  circuit  syn¬ 
thesis  problem  thus  reduces  to  one  of  synthesizing  terminal  eigenadmittances, 
such  as  that  expressed  in  (14) . 

The  synthesis  of  the  terminal  eigenadmittances  would  be  carried  out  in 
the  most  favorable  fashion  if  the  terms  in  the  summation  of  (15)  could  be 
dealt  with  a  conjugate  pole-pair  at  a  time.  Conjugate  pairs  of  modified 
pole  admittances  are  required  to  provide  the  conjugate  symmetry  associated 
with  PR  functions.  Synthesizing  these  terms  on  a  pole-pair  basis  limits  the 
constituent  admittances  which  are  dealt  with  to  second  order  in  s.  Unfor¬ 
tunately,  a  number  of  counterexamples  to  the  PR-ness  of  these  conjugate 
pole-pair  terms  are  described  in  reference  [11]. 

For  some  structures,  the  poles  associated  with  a  given  eigenvalue  com¬ 
prise  a  finite  collection.  An  analytically  tractable  example  of  a  structure 
manifesting  a  finite  collection  of  poles  associated  with  each  eigenvalue  is 
the  sphere  [3].  One  of  the  authors  of  the  present  paper  (DRW)  has  recently 
engaged  in  a  limited  study  of  the  eigenvalue /pole  associations  present  for 
the  straight  wire  scatterer.  The  evidence  uncovered  in  this  numerical  study 
of  the  eigenvalues  for  the  wire  Indicates  a  collection  of  poles  distributed 
along  an  arc  in  the  s-plane,  as  manifested  by  the  sphere.  On  the  other  hand, 
a  wire-loop  structure  evidences  an  infinite  number  of  poles  associated  with 
each  eigenvalue  [14].  The  sphere  and  wire  structures  may  be  differentiated 
from  the  loop  on  a  topological  basis:  Both  the  sphere  and  the  wire  are 
simply-connected  structures,  while  the  loop  is  a  multiply-connected  structure. 


17 


For  structures  where  the  pole  set  Is  finite  for  a  given  eigenvalue,  the 
summation  indicated  in  (14)  is  a  finite  summation.  This  resolves  the  trunca¬ 
tion  issue  in  the  summation  expressions  for  terminal  eigenadmittances.  In 
the  case  of  the  loop,  there  is  a  dominant  collection  of  poles  which  is  of 
finite  extent  (the  so-called  type  1  and  type  2  poles  [14]),  while  the  infinite 
collection  is  systematic  in  its  form  (the  so-called  type  3  poles).  The 
numerical  study  reported  in  reference  [11]  Indicates  that  the  sum  in  (14) 
converges  to  a  positive-real  value  for  a  reasonable  number  of  type  3  poles 
retained  in  the  summation  of  (14). 

Perhaps  more  important  than  the  truncation  issue  in  determining  the 
positive-realness  of  (14)  is  the  role  of  the  two  terms  Y^*  and  Y®nt.  There 
is  no  evidence  to  date  that  an  entire  function  contribution  ever  arises  in 
the  SEM  representation.  On  the  other  hand,  it  seems  likely  that  for  many 
geometries  branch  Integral  terms  do  appear  in  the  singularity  expansion  for 
reciprocal  eigenvalues.  When  the  eigenadmittances  are  all  summed  together, 
the  branch  integral  contributions  must  vanish  since  the  terminal  admittance 
is  known  to  be  meromorphic.  However,  the  branch  integrals  must  be  retained 
in  the  terminal  eigenadmittance  representation  in  order  to  insure  their 
PR-ness.  The  numerical  study  reported  In  reference  [11]  indicates  that  the 
terminal  eigenadmittance  derived  from  a  straight  wire  driven  at  its  one- 
quarter  point  manifests  a  non-PR  eigenadmittance  for  the  third  eigenset. 
Investigation  of  this  lack  of  PR-ness  indicates  that  the  departure  from  PR 
is  sufficiently  drastic  that  it  is  not  explainable  in  terms  of  numerical 
error  in  the  derived  SEM  parameters.  It  is  conceivable  that  the  failing  for 
this  particular  case  occurs  because  of  the  absence  of  the  branch  integral  in 
the  finite  pole  sum  taken  from  the  numerical  SEM  data.  The  issues  of  the 
occurrence  of  branch  Integrals  and  their  consequences  on  PR  properties  of 
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terminal  eigenadmittances  bear  further  investigation.  This  investigation  is 
substantially  complicated  by  the  fact  that  most  problems  which  are  treated 
by  means  of  the  Singularity  Expansion  Method  are  handled  on  a  numerical 
basis.  Indeed,  our  early  investigations  of  the  occurrence  of  branch  inte¬ 
grals  seem  to  indicate  that  their  presence  results  from  asymmetry  in  the 
structure.  Thus,  structures  which  are  highly  symmetrical,  such  as  the  sphere 
and  the  loop,  would  not  be  expected  to  yield  branch  integrals  in  thei'r 
analytically  obtained  eigenvalues. 

The  numerical  study  reported  in  [11]  indicates  that  for  both  the  straight 
wire  scatterer  and  for  the  loop,  the  single  dominant  resonant  pole— pair 
(i.e.,  the  pole-pair  located  nearest  the  jw  axis)  is  almost  PR  when  taken 
alone.  Furthermore,  this  single  pole  is  the  dominant  contributor  to  the 
summation  portion  of  the  terminal  eigenadmittance  represented  (14).  If  we 
choose  the  indices  +1  and  -1  to  denote  the  conjugate  constituents  of  the 
dominant  pole-pair  in  the  residue  expansion  for  the  reciprocal  eigenvalues 
and  for  the  terminal  eigenadmittances,  we  may  write 


Y  (s)  a  Yres(s) 
n  n 


<J_(r>s  .)  ;  a  > 
n  nl  S  S„ 

«  - a 

nl  .2 


r_i-  ♦  jj 

\-S  Snl  snlJ 


+  6 


-  -  *  2 

<J  (r,s  )  ;  a  >* 
n  nl  g  s 


nl 


b^r  *  -4-] 

is  —  s  ,  s  ,  J 
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(17) 


We  emphasize  that  the  presumption  of  this  approximation  is  an  empirical  one, 
based  on  observation  of  the  straight  wire  and  loop  results  alone.  The  occur¬ 
rence  of  a  significant  branch  integral  contribution  to  Y  is  certainly  one 

n 

\ 
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means  whereby  this  result  might  be  invalid.  For  the  numerical  examples 
tested  in  [91,  these  "almost  PR"  admittances  can  be  adjusted  to  be  approxi¬ 
mated  by  a  PR  admittance.  The  adjustment  required  is  negligible  in  mag¬ 
nitude  in  these  two  examples.  Thus,  we  are  able  to  write  the  approximation 

Y  (s)  *  G  +  Y^eS (s)  , 
n  n  n 

%  pes 

where  G  is  a  small,  positive  conductance  added  to  adjust  Y  to  be  PR.  If 
n  n 

this  approximation  is  possible  with  some  degree  of  generality,  it  constitutes 
a  particularly  robust  result.  Namely,  the  terminal  eigenadmittances  for  a 
structure  may  be  approximately  represented  in  terms  of  second-order  admit¬ 
tances  comprising  an  adjusted  pole-pair  admittance.  Such  a  second-order 
admittance  is  relatively  simple  to  realize. 

V.  APPLICABLE  CIRCUIT  SYNTHESIS  TECHNIQUES 

The  problem  of  translating  a  positive-real  rational  function,  such  as 
one  derives  for  a  pole-pair  or  for  an  eigenset  of  poles,  falls  in  the  realm 
of  classical  network  synthesis.  The  applicable  synthesis  methods  have 
existed  for  a  long  time  [15,16,17].  To  the  end  of  providing  a  realization 
of  equivalent  admittances  in  terms  of  RLC  networks,  we  present  a  brief  over¬ 
view  of  the  methods  available.  Although  there  exist  methods  which  include 
ideal  transformers  and/or  gyrators  in  the  realization,  we  avoid  these  realiza¬ 
tions  since  many  applications  preclude  the  use  of  these  devices  (e.g.,  due 
to  power-handling  considerations). 

The  classical  approach  to  the  synthesis  of  a  lossy  one-port  is  the  Brune 
method  [15,16].  The  Brune  method  realizes  the  sought-after  admittance  or 
impedance  in  terms  of  a  ladder  network.  This  method  was  significant  at  the 
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time  that  it  was  originally  formulated,  because  it  was  the  first  successful 
effort  to  realize  a  general  lossy  PR  one-port.  The  use  of  the  ladder  topology 
introduces  a  constraint  which  results  in  an  undesirable  element  in  the 
realization.  This  undesirable  element  appears  either  as  a  perfectly-coupled 
transformer  or  a  negative  inductance.  Neither  of  these  elements  is  realizable 
without  the  use  of  active  devices,  of  course.  Therefore,  the  Brune  procedure 
offers  little  benefic  for  present  purposes. 

The  Bott-Duffin  procedure  was  the  first  procedure  developed  which 
allowed  the  synthesis  of  a  general  PR  one-port  admittance  or  impedance  with 
a  network  comprising  only  resistors,  capacitors,  and  inductors.  The  Bott- 
Duffin  procedure  is  useful  for  formal  circuits  in  the  present  work  because 
of  its  generality  and  its  relative  amenability  to  automated  implementation. 

It  is  somewhat  doubtful,  however,  that  it  is  useful  for  the  synthesis  of 
implementable  circuits.  This  doubt  stems  from  the  fact  that  the  Bott-Duffin 
network  takes  the  form  of  a  balanced  bridge  for  all  frequencies  [16].  Indeed, 
there  Is  a  measure  of  arbitrariness  in  the  network  in  that  any  impedance 
which  one  chooses  may  be  inserted  in  the  so-called  detector  leg  of  the  bridge. 
That  the  network  constitutes  a  balanced  bridge  produces  a  high  measure  of 
sensitivity  of  the  network  response  to  errors  in  component  values  In  the 
network.  This  fact  is  dramatically  demonstrated  in  the  present  context  by  an 
example  given  in  reference  [11].  Balabanian  gives  some  transformations  of  the 
Bott-Duffin  topology  which  result  from  different  choices  of  the  arbitrary 
impedance  in  the  bridge  leg  of  the  network  [16].  It  is  conceivable  that  some 
or  all  of  these  transformations  can  provide  a  utilitarian  synthesis  by  sup¬ 
pressing  the  sensitivity  to  component  errors. 

The  Miyata  synthesis  procedure  [15,16]  offers  a  more  practical  synthesis 
method  than  either  the  Brune  or  the  Bott-Duffin  procedures.  The  circuits 
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which  result  from  this  procedure  are  always  RLC  and  would  seem  to  present 
less  sensitivity  to  component  errors.  The  implementation  of  this  method  is 
a  bit  more  complicated  than  the  implementation  of  either  of  the  previous 
two,  because  no  single  format  of  the  synthesis  is  applicable  to  a  general  PR 
function.  The  method  must  be  adapted  on  what  is  quite  nearly  an  ad  hoc 
basis  to  make  adjustments  in  the  particular  form  of  the  rational  function 
with  which  one  deals  [15]. 

The  fourth  alternative  available  as  a  synthesis  procedure  is  the 
Darlington  method  [17].  This  method  realizes  the  admittance  in  terms  of  a 
lossless  two-port  which  is  terminated  in  a  resistance.  Although  the 
Darlington  procedure  produces  networks  which  are  relatively  insensitive  to 
component  error,  it  has  the  drawback  that  it  may  or  may  not  result  in  a 
purely  RLC  network.  The  method  for  some  admittances  or  impedances  can  intro¬ 
duce  either  gyrators  or  ideal  transformers. 


VI.  SOURCE  SYNTHESIS  CONSIDERATIONS 

While  positive  real  considerations  constitute  the  principal  issue  in 
the  realizability  of  the  admittance  constituents  of  the  equivalent  circuit, 
we  must  consider  the  scheme  of  implementation  carefully  in  postulating  the 
means  whereby  the  associated  generator  elements  might  be  realized.  We  con¬ 
sider  two  generator  forms  here:  namely,  the  eigenvoltage  sources  given  by 
(11)  and  the  pole  voltage  sources  given  by  (16). 

Most  cases  of  practical  interest  for  broadband  circuit  synthesis  involve 
an  energy  collector  which  is  excited  by  what  we  shall  term  a  "factorable" 
incident  field.  That  is ,  a  field  which  can  be  written  in  the  form 

Einc(r,s)  -  F(r)  p(r,s)  f(s)  . 


(18) 
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Here,  F(r)  indicates  a  spatial  distribution  factor,  p(r,s)  is  a  propagation 
factor  which  links  space  and  time  dependencies,  and  f(a)  is  the  Laplace 
transform  of  the  time  history  of  the  wave.  For  example,  a  plane  wave  with 
a  time  history  which  is  a  Dirac  function,  6(t),  may  be  written 

Einc(r,s)  *  Eq  •  (I  -  nn)  e~sP*^c  i  . 

The  factor  1  is  the  transform  of  6(t),  Eq  is  a  constant  polarization  vector, 
and  p  is  a  unit  vector  in  the  direction  of  propagation.  If  the  plane  wave 
carries  a  time  history  f(t)  *  X  ^{f (s) } ,  then  the  representation  is 

flac<;,s)  -  E0  ■  (f  -  40  «.)  . 

A  highly  desirable  scheme  of  implementation  for  the  elgensource  form 
would  introduce  f(t)  into  only  a  single  port.  This  concept  is  understood 
more  clearly  with  a  reference  to  Figure  4.  A  voltage  v(t)  *  /it"  f(t)  is 
introduced  at  the  port  labeled  "INPUT”.  The  introduc  ion  of  f  (t)  at  a 
single  port  allows  the  use  of  the  network  with  any  input  waveform  which 
adheres  to  the  bandwidth  limitations  inherent  in  the  circuit  design.  The 
input  is  divided  into  N  identical  signals  f(t),  where  N  is  the  number  of 
constituent  eigensources.  These  sources  would  be  processed  by  transfer 
functions  Gn(s)  which  modify  the  spectrum  of  f(s)  to  produce  Vfl(s)  as  given 
by  (11).  This  transfer  function  is  identifiable  from  (11)  and  (18),  Viz., 

>  <J  (r,s)  ;  F(r)  p(r,s)> 

V  (s)  -  A  — -  f  (s) 

n  •  /*  x  a 

<Jn(r,s)  ;  ag> 

•  Gn(s)  f (s)  . 


(19) 
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In  ehis  implementation,  the  issue  becomes,  then,  the  realizability  of 

the  Gq(s).  Note  that  the  Gr  depend  on  the  source  configuration  that  produces 

the  incident  field — i.e.,  on  parameters  F  and  p.  If,  for  example,  in  the 

case  of  plane  wave  excitation,  the  angle  of  incidence  is  changed,  then  the 

transfer  functions  G  also  change. 

n 

If  the  circuit  may  be  limited  to  a  realization  of  the  energy  collecting 

properties  of  the  object  for  a  single  given  excitation  configuration,  then 

the  collection  of  transfer  functions  G  (s)  is  fixed.  It  is  conceivable  that 

n 

through  the  well  understood  approximation  and  synthesis  methods  for  two-port 
networks,  the  Gn(s)  might  be  realized,  at  least  in  an  approximate  fashion. 

Another  practical  possibility  is  the  realization  of  transfer  functions 
which  are  constant  amplitude  and  linear  phase  shift  with  respect  to  fre¬ 
quency.  That  is,  they  are  realizable  in  terms  of  an  attenuation  and  a  time 

delay.  Clearly,  such  a  realization  is  possible  only  if  the  G  (s)  in  (19) 

n 

are  independent  of  frequency,  at  least  in  an  approximate  sense.  We  shall 
explore  this  possibility  by  means  of  the  example  of  the  straight  wire  scat- 
terer. 

The  eigenmodes  of  the  straight  wire  scatterer  are  observed  to  be  essen¬ 
tially  independent  of  frequency  and  can  be  approximated  by  trigonometric 
functions.  For  the  coordinate  system  indicated  in  Figure  5,  we  write  the 
net  longitudinal  current  eigenmodes  on  the  wire  as 


In(z,s)  * 
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(20) 


Figure  5.  Incident  wave  geometry  for  the  straight  wire  sources 


26 


The  z  component  of  the  incident  plane  wave  is  expressible  on  the  wire  as 
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Thus, 
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Evaluation  of  the  integral  for  the  modes  approximated  by  (20)  yields 
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Recall  that  we  indicate  in  Section  XI  that  a  voltage  source  and  hence  the 

associated  G  is  not  needed  when  the  port  lies  at  a  null  of  the  mode.  Hence, 
n 

we  never  evaluate  (23)  in  a  situation  that  produces  a  zero  in  the  denomina¬ 
tor. 

It  is  clear  in  (23)  that  the  source  transfer  functions  are  significantly 
dependent  on  frequency  and  that  the  nature  of  this  frequency  dependence 
changes  with  the  incident  angle  0.  Only  in  the  broadside  incidence  case, 
0*0,  are  the  ^^(s)  frequency  independent.  This  indicates  in  a  concrete 
fashion  what  one  is  tempted  to  discern  from  (19):  The  frequency  dependence 
in  the  source  transfer  functions  stems  from  phase  delay  as  the  incident  field 
propagates  across  the  object.  This  dependence  exhibited  appears  physically 
realizable,  however:  the  polynomial  denominator  in  terms  of  a  resonant 
circuit  and  the  hyperbolic  functions  in  terms  of  time  delays.  The  time 
delays  are  recognizable  by  factoring  a  common  time  advance  and  shifting  the 
time  origin  to  suppress  this  advance. 
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Another  tempting  alternative  is  to  exploit  the  frequency  independent 
character  of  the  pole  voltage  sources  of  (16).  However,  the  admittances 
associated  therewith  are  realizable  on  a  conjugate  pair  basis.  Thus,  the 
two  pole  terms  must  be  combined.  We  adopt  the  notation  that  i  and  -i  con¬ 
stitute  the  indices  for  a  conjugate  pair  of  poles  for  a  given  eigenfunction. 
Thus,  the  conjugate  pair  admittance  is 


Yni(s)  "  Y-(s)  +  Y-  -*<s> 


ni 


n»_i 


and  the  pole  current  source  is 


IC?(s)  =  Y  (s)  V.  +  Yn  ,(s)  V  ,  . 
ni  ni  ni  n,-i  n,-i 


(24) 


Baum  [8]  points  out  that  only  for  the  case  that  V  .  is  real  such  that  V  .  » 

ni  ni 

V  t  does  a  frequency  independent  source  result.  He  also  points  out  that 

factoring  (24)  to  obtain  a  voltage  source  in  series  with  Yc^(s)  leads  to 

ni 

frequency  dependent  voltage  sources.  The  nature  of  the  frequency  dependence 

again  changes  with  the  incident  field  configuration  since  the  V  do. 

ni 

The  observations  and  comments  in  this  section  delve  into  only  the  more 
obvious  approaches  to  generator  synthesis.  The  approximation  and  synthesis 
of  the  transfer  function  form  In  (L9)  will,  possibly  prove  fruitful  for  a 
fixed  spatial  form  of  excitation.  The  simple  attenuator/delay  representation 
is  likely  to  be  useful  for  electrically-thin  objects*  when  the  incident  wave 
traverses  the  electrically-thin  direction. 


*Electrically-thin  in  terms  of  the  shortest  significant  wavelength  in  the 
excitation  spectrum. 
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VII.  CONCLUSIONS 

Several  Issues  are  discussed  in  this  work  relative  to  the  practical 
realizability  of  equivalent  circuits  of  energy  collecting  structures  where 
the  circuits  are  derived  from  the  singularity  expansion  of  the  scattering 
response  on  the  structure.  The  principal  issue  presented  is  whether  or  not 
component  impedance  or  admittance  functions  which  result  are  positive-real 
(and  hence  realizable)  functions.  A  general  result  is  that  the  terminal 
eigenadmittance  quantities  defined  herein  are  positive-real  functions  of 
frequency. 

We  introduce  the  concept  of  approximating  terminal  eigenadmittances  by 
the  pole  admittance  of  the  most  highly  resonant  pole  in  each  eigenset.  This 
approach  is  supported  as  to  its  practicability  in  reference  [11]  through  the 
examples  of  the  straight  wire  scatterer  and  the  wire  loop.  It  is  emphasized 
that  this  approximation  is  ad  hoc  in  character  and  must  be  tested  for  any 
specific  structure  to  which  it  is  applied.  The  cases  cited  produced  admit¬ 
tances  which  were  "almost  positive-real" — which  could  be  shifted  negligibly 
to  be  made  positive-real.  The  quality  of  the  approximation  and  the  near 
PR-ness  of  the  resulting  admittance  is  conceivably  a  consequence  of  the  high 
Q  character  of  the  structures  considered.  For  a  low  Q  structure  such  as  a 
spherical  scatterer,  it  is  not  clear  whether  or  not  the  admittance  contribu¬ 
tion  of  the  dominant  resonance  pole  in  the  eigenset  will  suffice  or  prove 
realizable. 

Some  of  the  practical  considerations  in  implementing  the  source  portions 
of  the  equivalent  circuits  are  discussed.  In  particular,  the  desirability  of 
formulating  aspect  and/or  frequency  independent  transfer  functions  from  the 
incident  waveform's  time  history  to  the  individual  sources  is  suggested. 
Indications  are  chat  each  of  these  features  results  only  when  geometric 


degeneracies  are  present  and  are  not  true  in  general.  If  the  frequency  and 
aspect  dependent  character  of  these  transfer  functions  Is  borne  out,  the 
upshot  is  that  the  transfer  functions  must  be  realized  through  the  approxima¬ 
tion  and  realization  procedures  for  two-port  circuit  synthesis.  This  realiza¬ 
tion  will  be  specific  to  a  given  angle  of  incidence.  Although  a  more  general 
result  is  desirable,  in  many  applications,  such  as  the  modeling  of  energy 
coupling  into  a  scatterer  which  functions  as  a  shield  as  well,  the  single- 
aspect-angle  circuit  can  provide  a  useful  worst  case  model. 

It  is  suggested  herein  that  for  some  structures—* in  particular  those 
devoid  of  appreciable  symmetry — branch  integrals  may  constitute  a  portion  of 
the  terminal  eigenadmittance  representation  along  with  the  usual  SEM  pole 
terms.  Though  these  branch  integrals  Ho  not  appear  in  the  ultimate  current 
representation,  they  can  well  compromise  the  realizability  of  the  eigen¬ 
admittance  segments  into  which  the  circuit  realization  is  parceled.  The 
evidence  available  regarding  these  branch  integrals  is  sketchy,  at  best. 
Further  work  is  warranted  in  this  area,  too. 

In  summary,  the  practical  realizability  of  circuits  which  model  the 
energy  collecting  properties  of  equivalent  circuits  over  a  broad  range  of 
frequencies  appears  possible.  In  candor,  the  procedures  and  understanding 
of  realizability  issues  for  a  general  structure  must  be  declared  primitive, 
at  present. 

VIII.  ACKNOWLEDGMENTS 

The  authors  wish  to  acknowledge  the  many  fruitful  discussions  on  this 
subject  with  Dr.  Carl  Baum  of  the  Air  Force  Weapons  Laboratory  and  with 
Mr.  Gary  Streable,  now  with  Bell  Laboratories.  Some  of  the  material  in 
Section  V  stems  from  a  discussion  with  Prof.  Norman  Balabanian  of  Syracuse 


University.  Ms.  Gail  McAlister  provided  the  capable  typing  and  editing  of 
the  manuscript. 


31 


Appendix 

Consider  a  conducting  scatterer  S  with  unit  surface  normal  n  illuminated 

p  f*  _  a  a 

by  an  incident  tangential  electric  field  E  (r,t).  A  surface  current  J(r,t) 
is  induced  on  S,  producing  a  scattered  electric  field  Esc(r,t)  which  satis¬ 
fies  the  boundary  condition 

ElllC(r, C)  +  '  0  ,  rcS.  (A.1) 

By  the  equivalence  theorem,  the  scatterer  may  be  replaced  by  a  surface  cur¬ 
rent  J(r,t)  on  S  radiating  in  free  space.  This  current,  when  radiating  in 
the  absence  of  the  sources  of  the  incident  field,  produces  the  scattered 
field  Esc  exterior  to  S. 

In  the  absence  of  other  sources,  the  total  energy  U  radiated  by  the  cur¬ 
rent  distribution  J  is  positive-semidefinite; 

t 

U(t)  -  f  W(t)  dt  >  0  (A. 2) 

—09 

where  U  is  the  energy  and  W(t)  is  the  total  power  radiated  by  J  at  time  t  and 
is  computed  as  follows: 

W(t)  -  -<I8ctf,t)  ;  J(r,t)> 

-  <Einc(?,t)  ;  J(r,t)>  .  (A. 3) 

•i  «Vi  *4  ftp 

Although  the  current  J(r,t)  and  the  incident  field  E  (r,t)  are  related  by 
an  integral  equation  (1) ,  we  do  not  make  use  of  this  fact  until  later. 
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Next,  suppose  the  incident  electric  field  and  induced  current  have  the 


form 


J(r,t)  *  J(r,s)  est  +  J*(r,s)  es  C 


st. 


2  Re{J(r,s)  e  >  , 


(A. 4a) 


and 


Einc(?,t) 


* 

Einc(;,s)  eSt  +  Einc*(?,s)  eS  C 
2  Re{f±nc (r,s)  eSt}  , 


(A. 4b) 


■■ 

where  s  •  a  +  ja>  is  a  fixed  complex  number  with  a  >  0,  E  (r,s)  is  an  arbi¬ 
trary  complex  vector  function  of  position,  and  J(r,s)  is  the  spatial  dis¬ 
tribution  of  the  resulting  complex  current  response.  The  excitation  is 
assumed  to  begin  at  t  B  at  which  time  there  is  no  initial  energy  in  the 
system.  Since  eat  *  0  for  t  *»  -*•,  both  Einc  and  J  are  zero  at  t  »  -*•  and 
there  is  no  transient  term.  That  is,  the  forced  response  (A. 4a)  is  the  total 
response. 

The  power  radiated  W(t)  is  now 


W(t) 


+ 


<3(i,s)  ;  linc(?,s)>  e2st 

A 

,  sine*,-  .  2s  t 
<J  (r,s)  ;  E  (r,s)>  e 


+  [<Einc(r,s)  ;  J  (r,s)>  +  <Einc*(r,s)  ;  J(r,s)>]  e2at 
-  2  Re{<finc(r,s)  ;  3(r,s)>  e2st 
+  <l(?,s)  ;  J*(J,s)>  e2at}  . 
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Hence,  the  total  energy  radiated,  from  (A. 2),  is 

U(t)  -  Re{<linc(r,s)  ;  J(r,s)>  e2st/s 

+  <E*nc(r,s)  ;  J*(r,s)>  e2ot/0>  . 

If  we  write  the  first  of  the  above  products  in  polar  form, 

<Einc(r,s)  ;  J(r,s)>/s  =*  |<Einc(r,s)  ;  J(r,s)>/sj  , 
then  U(t)  can  be  written  as 

U(t)  =  e2ot  ( |  <EJnc(r,s)  ;  J(r,s)>/s|  cos(2mt  +  4>) 

+  Re{<Einc(r,s)  ;  J*(r,s)>/o}]  >  0  . 

Me  consider  the  cases  w  =  0  and  <o  ^  0  separately. 

CASE  I ,  w  *  0 . 

U(t)  =  2  <Re{EitlC(r,s) }  ;  Re{J(r,s)}>  e2ot/o  , 

whence 

<Re{linC(r,s) }  ;  Re{jU,s)>>  >  0  . 

s^nc  »  “* in c  • 

Since  E  (r,s)  is  arbitrary,  we  could  replace  it  with  jE  (r,s 

corresponding  current  response  would  be  jJ(r,s).  Hence,  from  (A 

follows  that 

<Re{jEinc(r,s) }  ;  Re{jJ(r,s)}>  «  <Im{Einc(r,s) }  ;  Im{J(r,s)}: 
Equations  (A. 7)  and  (A. 8)  imply  that 


(A.  5) 


(A.  6) 


(A.  7) 

and  the 
7)  it 

>  0  .  (A.  8) 


Re(<E^nc(r,s)  ;  J  (r,s)>>  >  0  . 


(A.  9) 
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CASE  II,  ui  i  0. 

Since  mln{cos(o)t  +  ♦)}  ■  -1,  we  conclude  from  (A.  6)  that 

Re{<linc(r,s)  ;  3*(r,s)>}/o  >  |<finc(?,s)  ;  3(^,s)>/s|  >  0  .  (A.10) 

Thus,  we  again  obtain 

Re{<Einc(r,s)  ;  J*(r,s)>}  >  0  ,  a  >  0  (A.ll) 

as  a  necessary  condition  on  the  current  on  a  passive  scatterer. 

Using  (A.ll),  we  may  derive  a  condition  similar  to  the  positive-real 
condition  on  driving-point  immittances  in  circuit  theory.  The  condition 
applies  to  the  dyadic  kernel  of  (1).  Recognizing  that  the  quantities 

mm  mm  « 

E  (r,s)  and  J(r,s)  in  the  foregoing  are  the  transformed  field  quantities 
in  (1) ,  we  write 

Re{<J*(r,s)  ;  ?(r,r',s)  ;  J(r,s)> 

=  <J*(r,s)  ;  Re{r(r,r' ,s)}  ;  J(r,s)>  >  0  .  (A.12) 

This  result  is  obtained  with  the  benefit  of  the  dyadic  reciprocity  condition 

?+(r,r’,s)  -  ?*(r,  rf  ,s)  =*  f*(r?,r,s)  . 

Thus,  Re{r(r,r’ ,s) }  is  positive-semidefinite  for  a  >  0.  The  condition  (A.12) 
is  the  operator  counterpart  to  one  of  the  so-called  positive-real  conditions 
for  lmmlttance  matrices  of  multiport  networks.  We  may  establish  counterparts 
to  the  remaining  two  conditions  defining  PF-ness  for  multiports  by  examining 
the  dyadic  kernel  ?.  The  conditions  may  be  stated  as  follows: 

a  _  _ 

(i)  r(r,r',s)  is  analytic  for  Re{s}  >  0; 

(li)  r(r,r’,s)  is  real  for  real  positive  s; 


with  (A. 12)  restated  as 


(iii)  Re{?(r,r’ ,s)}  >  0  ,  Re{s)  >  0. 

The  immittance  matrix  counterparts  to  these  conditions  are  necessary  and 

sufficient  that  rational  polynomial  immittance  matrices  completely  represent 

a  network  comprising  passive  elements.  While  we  presently  cannot  establish 

Si  —  — 

sufficiency,  the  three  conditions  above  are  necessary  for  r(r,r',s)  to 
represent  a  passive  scattering  process. 

Not  only  the  dyadic  kernel,  but  also  any  input  immittance  quantity  is 
positive-real  when  the  object  is  excited  as  an  antenna.  To  see  this,  we 
begin  with  (A. 11)  with  the  excitation  specialized  to  a  field  impressed 


across  a  gap  surface 


namely. 


C  a  V(s) /A  ,  r  e  S 

Einc(r,s)  -  |  ^  8 

Lo  ,  ^Sg. 


(A. 13) 


Then  (A. 11)  becomes 

Re{V(s)  <a  ;  J  (r,s)>  /A}  =  Re{V(s)  I  (s)} 

R  Sg 

=  Re { V ( s )  Y*(s)  V* (s) } 

=  j V(s) | 2  Re{Y(s) }  >  0  ,  a  >  0  .  (A. 14) 


This  implies  that 

Re{Y(s)}  >  0  ,  rr  >  0  .  (A. 15) 

The  remaining  two  positive-real  conditions  follow  from  observation.  Thus  the 
input  admittance  from  the  short  circuit  boundary  value  problem  is  positive- 
real.  Similarly,  one  may  show  that  the  input  impedance  obtained  from  the 
open  circuit  boundary  value  problem  [8J  is  also  positive-real. 


36 


a 

Next,  suppose  that  the  dyadic  kernel  r  has  an  eigenvalue  Xq  and  corre¬ 
sponding  eigenvector  JQ  satisfying  (5).  If  we  choose  as  an  excitation 

finc(r,s)  -  Xq(s)  3nU,s)  , 

the  response  current  is  Jn(r,s)  and  (A. 11)  becomes 

Re{Xn(s)  <JQ(r,s)  ;  J*(r,s)>} 

■  Re{Xn(s)>  jjJn(r,s)||^  >  0  ,  a  >  0  . 

This  implies  that 

Re{X  (s)}  >  0  ,  a  >  0  .  (A. 16) 

n 

By  observation,  the  remaining  PR  conditions  may  be  established  for  the  eigen¬ 
values  with  the  exception  of  the  case  of  a  X^Cs)  with  a  branch  point  in  the 
right  half  plane.  We  may  admit  even  such  an  eigenvalue  to  an  "extended"  PR 
class  where  the  "extended"  implies  "analytic  except  for  branch  points"  in  the 
right  half  plane. 

Thus,  eigenvalues  (eigen impedances)  are  positive-real.  Since  recipro¬ 
cals  of  positive-real  quantities  are  also  positive-real,  then  the  reciprocal 
eigenvalues  (eigenadmittances)  must  also  be  positive-real; 

Re  ^  7  y  >  0  ,  a  >  0  .  (A. 17) 

Vs' 

For  scalar  positive-real  quantities,  such  as  input  immittances  or  elgen- 
immittances,  it  is  well-known  that  the  positive-real  conditions  are  also 
equivalent  to  the  following  conditions  [18]: 


(a)  f(s)  has  no  poles  or  zeros  for  o  >  0; 

(b)  poles  of  £(s)  on  the  Imaginary  axis  must  be  simple  and  their 
residues  must  be  real  and  positive;  and 

(c)  Re{f(jw)}  2:  0  ,  we  (0,»), 

where  f(s)  is  an  arbitrary  PR  function.  The  last  of  these  conditions  is 
particularly  important  in  testing  PR-ness  because  it  requires  that  f(s)  be 
examined  only  on  the  jw  axis. 
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